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Abstract 

Universal bounds on the electrical and elastic response of two-phase (and mul- 
tiphase) ellipsoidal or parallelopipedic bodies have been obtained by Nemat-Nasser 
and Hori. Here we show how their bounds can be improved and extended to bodies 
of arbitrary shape. Although our analysis is for two-phase bodies with isotropic 
phases it can easily be extended to multiphase bodies with anisotropic constituents. 
Our two-phase bounds can be used in an inverse fashion to bound the volume frac- 
tions occupied by the phases, and for electrical conductivity reduce to those of 
Capdeboscq and Vogelius when the volume fraction is asymptotically small. Other 
volume fraction bounds derived here utilize information obtained from thermal, 
magnetic, dielectric or elastic responses. One bound on the volume fraction can be 
obtained by simply immersing the body in a water filled cylinder with a piston at 
one end and measuring the change in water pressure when the piston is displaced 
by a known small amount. This bound may be particularly effective for estimating 
the volume of cavities in a body. We also obtain new bounds utilizing just one pair 
of (voltage, flux) electrical measurements at the boundary of the body. 

Keywords: size estimation, universal bounds, volume fraction bounds. 

1 Introduction 



Berryman and Kohn (1990) found that classical variational principles could be used to 
obtain information about the conductivity inside a body from electrical measurements 
on the exterior. In this paper our main focus is on using classical variational principles 
and known bounds on the response of periodic composites to bound the volume fraction 
of one phase in a two-phase body Q from measurements on the exterior of the body. Of 
course if one knows the mass densities of the two phases, the easiest way to do this is 
just to weigh the body. However this may not always be practical, or the densities of the 
two phases may be very close. 

Two types of boundary conditions are most natural: what we call special Dirichlet 
conditions where affine Dirichlet conditions are imposed on the boundary of Q (which 
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would render the field inside Q uniform if the body were homogeneous) or what we 
call special Neumann conditions where Neumann conditions are imposed that would 
render the field inside Q uniform if the body were homogeneous. Bounds on the elec- 
trical and elastic response of the body to these special boundary conditions were ob- 
tained by Nemat-Nasser and Hori (119931 fT995ll . and were extended to piezoelectricity by 



Hori and Nemat-Nasser (1998). They called these bounds universal because they did not 
depend on any assumption about the microgeometry in the body. They obtained both 
elementary universal bounds based on the classical variational principles, and reviewed 
below in section 3, and universal bounds based on the Hashin-Shtrikman variational 
principles (Hashin and Shtrikman I1962[ I1963p . The latter bounds were obtained under 
the assumption that Q is either an ellipsoid or a parallelopiped, but we will see here 
that they can easily be improved and generalized to bodies Q of arbitrary shape. The 
key is to consider an assemblage of copies of Q packed to fill all space, and then to 
use the bounds of Huet (1990) which relate the effective tensor of this composite to the 
responses of Q under the special boundary conditions. Then existing bounds on the effec- 
tive tensor [as surveyed in the books of Nemat-Nasser and Hori (1993), Cherkaev (2000), 



Allaire (2002| , [Milton (2002| ), and [Tartar (2009D ] can be directly applied to bound the 
responses of fl under special boundary conditions (see sections 5,6,7, and 8). Since 
these bounds involve the volume fractions of the phases (and the moduli of the phases), 
they can be used in an inverse fashion to bound the volume fraction. As shown by 



Kang, Kim, and Milton (2011) the volume fraction bounds thus obtained for electrical 



conductivity generalize those obtained by Capdeboscq and Vogelius (I2003[ I2004p for the 
important case when the volume fraction is asymptotically small. 

Given the close connection between bounds on effective tensors and bounds on the 
responses of Q under special boundary condition, a natural question to ask is whether 
methods that have been used to derive bounds on effective tensors could be directly 
used to derive bounds on the response of Q under more general boundary conditions. 
One such method is the Hashin-Strikman (I1962^ 119631) variational method and this lead 
Nemat-Nasser and Hori to their bounds for ellipsoidal or parallelopipedic Q. Another par- 
ticularly successful method is the translation method ([Tartar 1979[ Lurie and Cherkaev 
\WW2[ IMurat and Tartar 1985| ITartar 1985| IMilton 1990|l and indeed as shown in 

a companion paper ( Kang, Kim, and Milton 2011 ) this method yields upper and lower 
bounds on the volume fraction in a two-phase body with general boundary conditions for 
two-dimensional conductivity without making any assumption on the shape of Q. For 
special boundary conditions the bounds thus derived reduce to the ones derived here. 

We also provide (in section 4) some new conductivity bounds which just involve the 
results of just one (flux, voltage) pair measured at the boundary of Q, and which im- 
prove upon the elementary bounds of Nemat-Nasser and Hori (1993). Again these new 
bounds can be used in an inverse fashion to bound the volume fraction. Other volume 
fraction bounds using one measurement were derived by Kang, Seo, and Sheen (1997) 



Ikehata (1998), Alessandrini and Rosset (1998), Alessandrini, Rosset, and Seo (2000), and 



Alessandrini, Morassi, and Rosset (2002). These other bounds involve constants which 



are not easy to determine, making it difficult to make a general comparison with our new 
bounds. 

The various bounds on the volume fraction we have derived are too numerous to 
summarize in this introduction. However we want to draw attention to the bounds (I5.15P 
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and f l5.24p which are the natural extension of the famous Hashin and Shtrikman (1962) 
conductivity bounds to this problem. Also of particular note is the bound (I8.26p . which is 
one natural generalization of the bulk modulus bounds of Hashin and Shtrikman (1963) 
and Hill (1963), and implies that a bound on the volume fraction can be obtained by 
simply immersing the body in a water filled cylinder with a piston at one end and 
measuring the change in water pressure when the piston is displaced by a known small 
amount. 



2 The conductivity response tensors with special Dirich- 
let and special Neumann boundary conditions 

In electrical impedance tomography in a body Q containing two isotropic components 
with (positive, scalar) conductivities ui and a2 the potential V satisfies 

V ■ aVV = 0, where a(x) = x(x)(Ti + (1 - x(x))fT2, (2.1) 

and x(x) is the indicator function of component 1, taking the value 1 in component and 
in component 2. Equivalently, in terms of the current field j(x) and electric field e(x) 
we have 

V-j = 0, j = ae, e=-VK (2.2) 

Let us assume the components have been labeled so that ai > o"2. We are given a set of 
Cauchy data, i.e. measurements of pairs (Vo,g), where Vo(x) and g(x) are the boundary 
values of the voltage ^(x) and and flux g(x) = — n ■ j(x) at the boundary dQ of Q, in 
which n(x) is the outwards normal to the boundary. From this boundary information 
we can immediately determine, using integration by parts, volume averages such as 

l"l Jdn l"l Jan 



an 



(j) = ^ / -xg, (2.3) 



\n\ 



on 



where the angular brackets denote the volume average, i.e. 



for any quantity 5'(x). From such averages our objective is to bound the volume fraction 
fi = (x) of component 1 (and hence also the volume fraction /2 = 1 — /i of component 
2). 

To obtain good estimates of the volume fraction fi it makes physical sense to use 
measurements where the fields e(x) and j(x) probe well into the interior of Q. In this 
connection two sets of measurements are most natural. We could apply special Dirichlet 
boundary conditions 

Vo = -eo ■ X, (2.5) 
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and measure jo = (j). Here, according to fl2.3p . Gq equals (e). Since jo is linearly related 
to eo we can write 

jo = T^eo, (2.6) 

which defines the conductivity tensor cr^ [D for Dirichlet). To determine cr^ in dimen- 
sion d = 2, 3 it of course suffices to measure jo for d linearly independent values of eo. 
Alternatively we could apply the special Neumann boundary conditions 

g=jo-n, (2.7) 

and measure eo = (e). Again according to (12. 3p . jo = (j). The linear relation between eo 
and jo, 

eo = (^'^)-'jo (2.8) 

defines the resistivity tensor (cr^)~^ and hence the conductivity tensor cr^ {N for Neu- 
mann): we will see later that (cr^)~^ is invertible. To determine cr^ it suffices to measure 
eo for d linearly independent values of jo- With either of these two sorts of boundary 
conditions (but not in general) Hill (1963) has shown that 

(e-j) = (e)-(j), (2.9) 

as follows by substituting (12. 5p or (12.70 in the first of equations (12. 3p . Using this rela- 
tionship, and its obvious generalizations, it is easy to check that both cr^ and cr^ are 
self-adjoint. Thus if e'(x) and j'(x) denote the electric and current fields associated with 
the boundary conditions (12.51) . with eo replaced by eg, while keeping the same conduc- 
tivity cr(x) then 

e[, ■ ^T^eo = (e' ■ j) = (eVe) = (e ■ j') = eo • a^'e'^, (2.10) 
which implies cr^ is self-adjoint. By similar argument cr^ is self-adjoint. 



3 Known elementary bounds 



This section reviews the elementary bounds on cr^ and cr^ obtained by Nemat-Nasser and Hori (1993) 
and by Willis in a 1989 private communication to Nemat-Nasser and Hori. Their im- 
plications for bounding the volume fraction will be studied. We will make use of two 
classical variational principles: the Dirichlet variational principle that 



min / e • ae = / — Vog, (3.1) 

e(x) = — VVl(x) Jn Jan 

V{x.) = Voix.) on dn 

which is attained when V_{x.) = V{x.), and the Neumann variational principle that 

mill / j ■ a-'} = I Voq, (3.2) 

j(x) Jo. Jan 

V-[(x) = 

n • j(x) = — g(x) on dfl 
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which is attained when j(x) = j(x). With the special Dirichlet boundary conditions (12 .Sp 
the Dirichlet variational principle implies 

min (e ■ ae) = en • cr^en. (3.3) 
e(x) = -VZ(x) " 
V_{'x) = —Go ■ X on dVl 

Taking a trial potential V_ = — Gq • x produces the elementary upper bound on Gq • cr^eo 

eo ■ o-^eo < (o-)eo ■ eo, (3.4) 



given by Nemat-Nasser and Hori (1993). To obtain a lower bound observe, following a 
standard argument, that the left hand side of (13.31) is surely decreased if we take the 
minimum over a larger class of trial fields. Since the constraints on e(x) imply (e) = eg 
let us replace them by this weaker constraint to obtain the inequality 

eo • cr^eo > min (e-ae), (3.5) 
e(x) 

(e) = eo 

where the minimum is now over fields e which are not necessarily curl-free. Using La- 
grange multipliers one finds that the minimum is attained when e = a~^{a^^)^^eo and 
so we obtain the lower bound 

eo ■ cr^eo > (o-"^)"^eo ■ eo (3.6) 



of Nemat-Nasser and Hori (1993). Taken together, (13. 4p and (13. 6p imply the lower and 
upper bounds 

/ej^ - a,) /(a, - a,) < h < - -p^) - (3.7) 

V eo-eo J V eo-cr^eo/ 

on the volume fraction /i. These bounds give useful information even if we only know 
jo = CT^eo for only one value of eo, i.e. if we only take one measurement. These bounds 
(13. 7p are sharp in the sense that the lower bound is approached artitrarily closely if VL is 
filled with a periodic laminate of components 1 and 2, oriented with the normal to the 
layers orthogonal to eo and we let the period length go to zero, while the upper bound 
is approached artitrarily closely for the same geometry, but oriented with the normal to 
the layers parallel to eo. If the full tensor cr^ is known, from c? = 2, 3 measurements of 
pairs (eo, jo) then we can take the intersection of the bounds (13. 7p as eo is varied, and so 
obtain 

(A^ - (72)/(ai - (72) < /i < (1/^2 - l/A?)/((72-i - af^), (3.8) 

where A;^ and A^ are the maximum and minimum eigenvalues of cr^. However we will 
see in the next section that an additional and typically sharper upper bound on /i can 
be obtained. 

With the special Neumann boundary conditions (12. 7p the variational principle (13. 2 p 
implies 

min (j-a"'j)=jo-(0"'jo. (3.9) 

j(x) 
V-j(x) = 
n • j(x) = n ■ jo on 9f2 
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By taking a constant trial field j(x) = jo or alternatively by taking the minimum over 
the larger class of trial fields satisfying only (j) = jo we obtain the bounds 



^jo ■ jo < jc 



.N\-l 



of Nemat-Nasser and Hori (1993) which imply 



Jo ■ Jo 



jo • (^^)-^jo 



0-2 /(O"! - 0"2) < fl< CT; 



jo < {(^ ^)jo ■ jo 



-1 jo-('T^)-^jo 



Jo ■ Jo 



(3.10) 



(3.11) 



These bounds are applicable even if we only know Gq 



jo for only one value of 



jo- For comparison, with these special Neumann boundary conditions (12.71) . the bounds 



in Theorem 3.1 of Kang, Seo, and Sheen (1997) coupled with the improvement in propo- 
sition of Ikehata (1998), with cti = /c > 1, (T2 = 1 and jo • jo = 1; imply 



l-jo-(0-'jo)</i< 



k-l 



k 



k 



;i - jc 



(-^)-^jo). 



(3.12) 



In this case it is easy to check that the upper bounds in (13. lip and (I3.12p coincide while 
the lower bound in (I3.1ip is tighter. The bounds (13. lip are each approached arbitrarily 
closely if Q is filled with a periodic laminate of components 1 and 2, oriented with jo 
either parallel or orthogonal to the layers and we let the period length go to zero. 
In summary, (l3.4p .f lXB]) and (l3.1Up imply the matrix inequalities 



{(T-')-^I <ct'' < (a)I, {(J-T^I <a^ < {a)I 



(3.13) 



of Nemat-Nasser and Hori (1993). 

For artitrary boundary conditions, i.e. for any e and j satisfying (12. 2p within Q, we 
have the bounds 

(3.14) 



(e-j) > eo-cr^eo, (e-j)>jo(cr 



D\-l 



JO; 



where eo = (e) and jo = (j), due to Willis in a 1989 private communication to Nemat- 



Nasser and Hori, and presented by Nemat-Nasser and Hori (1993). In conjunction with 
(I3.13P they imply the volume fraction bounds. 



Jo ■ Jo 



(e-j) 



- - 0-2 /((Ti - (T2) < /l < {(^2^ 



eo ■ eo 
(e-j) 



(3.15) 



4 New bounds with one measurement 

If we have measurements of (e ■ j) and both vectors eo and jo for arbitrary boundary 
conditions then the bounds (13.141) and (I3.15P can be improved. The classical variational 
principle (13. ip implies 



min (e ■ ere), = (e ■ i) 

e(x) = -VHx) 

y(x) = 1/0 (x) on dn 

(f^e) = jo 



(4.1) 
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where we have chosen to add the constraint that (ere) = jo since we know that without 
this constraint the minimizer e = e satisfies (ae) = jo- We surely obtain something lower 
if take the minimum over the larger class of fields satisfying only (e) = Bq and (ae) = jg. 
Thus we obtain the inequality 

min (e ■ ae) < (e ■ j). (4.2) 
e(x) 

(e) = eo 
(o-e) = jo 

By introducing two vector valued Lagrange multipliers associated with the two vector 
valued constraints we find that the minimum is attained when 

e(x) = eo + {{a-') - a-\^)){{a){a-') - l)-^(jo - (a)eo). (4.3) 

Substituting this back in (14.21) gives the bound 

{{a) {a-') - l)((e ■ j) - eo ■ jo) > (jo - (^)eo) ■ ((^-')jo - eo). (4.4) 

If, with general boundary conditions, we are interested in bounding the volume fraction 
/i given measured values of (e ■ j), eo and jo then the difference between the left hand 
side and right hand side of (14. 4p is a quadratic in /i whose two roots give upper and 
lower bounds on /i. (Unless the roots happen to be complex, in which case there is no 
configuration of the two phases within fl which produce the measured (e • j), eo and jo, 
indicating the presence of other phases or indicating an error in measurements.) 

In the particular cases of either special Dirichlet or special Neumann boundary con- 
ditions, (12.51) or (12.71) . the left hand side of (14. 4p vanishes (see (12. 9p ) and we obtain the 
reduced bounds 

0>(jo-(a)eo)-((a-i)jo-eo), (4.5) 

which are in fact implied by the matrix inequalities (I3.13p . This bound (14. 5 p is optimal. 
For any given fixed eo, and fixed volume fraction /i, the vector jo has an endpoint 
which is constrained by (14. 5 p to lie within a sphere (disk in two dimensions) centered at 
(('^) + ('^~^)~^)eo/2. When Q is filled with a periodic laminate of the two phases with 
interfaces orthogonal to some unit vector m, and we let the period length go to zero, then 
the endpoint of the vector jo covers the entire surface of this sphere (disk) as m ranges 
over all unit vectors. These bounds are the analogs, for arbitrary bodies Q, of bounds on 



possible (eo, jo) pairs for composites derived by Raitum (1983) and Tartar (1995). If we 
are given eo and jo and want to bound /i then we should find the range of /i where the 
sphere (or disk) contains the endpoint of the vector jo. The endpoints of this range are 
the roots of the right hand side of (14. 5 p which is a quadratic function of /i. 

Knowledge of eo and jo is equivalent to knowledge of (e ■ v) and (v ■ j) for all constant 
fields V. A more general alternative is to use the information about 

fli = (e ■ j») = 7^ / -Vo{ii-n), 



'an 

6fe = (Vyfe-j)= / -Vk{yn), (4.6) 
l"l Ian 
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for a given set of "comparison flux fields" ji(x) satisfying V ■ jj = 0, i = 1, 2, . . . n and 
"comparison potentials" Vfc(x), k = 1,2, . . .m. Suppose, for example, that we have just 
one comparison flux field ji. We have the variational principle 

, , ^ (e-ae) = (e-j), (4.7) 

e(x) = -VV:(x) 
V:(x) = 1/o(x) on dQ 
(e ■ ji) = ai 

where we have chosen to add the constraint that (e- ji) = oi since we know that without 
this constraint the minimizer e = e satisfies (e ■ ji) = ai. This implies the inequality 

(e-j)> min (e-cre). (4.8) 
e(x) 

(e • ji) = ai 

By introducing a Lagrange multiplier associated with the constraint (e ■ ji) = ai we see 
the minimum occurs when 

e = ai(T-iji/(ji -a-iji), (4.9) 

giving the inequality 

{e-j)>al/{h-cT-%). (4.10) 

This inequality gives information about cr(x) through (ji -cr^^ji). If we only want bounds 
which involve the volume fraction we should choose ji(x) with 

|ji(x)| = l forallxGfi. (4.11) 

There are many divergence free fields ji(x) which satisfy this constraint. For example in 
two dimensions we can take 

= {d(f)/dx2, -90/(9x1 ), with |V0(x)| = 1 for all x G 1^. (4.12) 

Thus 0(x) satisfies an Eikonal equation, and we could take 0(x) to be the shortest 
distance between x and a curve outside Q. Once (14. lip is satisfied (14.101) implies the 
volume fraction bound 

,-1 1 \ /f^-i 



/i < [^^2' - /i^2' - ^i)- (4.13) 

In the special case when ji = eo/|eo| this reduces to the upper bound on /i given by 
(EISD. 

An important question is whether this new bound is sharp, and if so for what cr(x)? 
The new bound will be sharp when e = e where e is given by (14. 9p . In that case 

j(x) = aiji/(ji-a-iji) (4.14) 

has zero divergence because it is proportional to ji. Let us impose the Neumann boundary 
condition 

j(x) • n = ji ■ n for all x e dil, (4.15) 
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Figure 1: A schematic of the type of layered microstructure achieving the volume fraction 
bound fl4.13p . where the black regions denote one phase, and the white regions the other 
phase. The layer widths should be much finer than the size of Q. 

and look for a o"(x) so j(x) = ji(x) and e(x) = (T~^ji(x) is curl-free. Now as schematically 
represented by figure [H choose cr(x) to correspond to a finely layered composite with 
layers orthogonal to the streamlines of ji(x), and with phase 1 occupying a local volume 
fraction p(x). This composite will support a current field j(x) = ji(x) and an electric 
field e(x) = o-^^ji(x) provided 

V X eo = 0, eo = [a^' - pix){a^' - a^')]iii^). (4.16) 

Here eo(x) is the weak limit (local volume average) of e(x) as the layer spacing goes 
to zero. In two dimensions, given ji(x) we could look for solutions for p(x) such that 
(14.161) is satisfied and < p(x) < 1 in f2. We expect such solutions to exist for a wide 
class of fields ji(x). This example shows that non-constant "comparison flux fields" can 
lead to sharp bounds on the volume fraction. In three dimensions we only expect to find 
a solution of the vector equation (14.161) for the scalar field p(x) if ji(x) satisfies some 
additional conditions. 

5 Relationship to bounding effective tensors of com- 
posites 

Consider a periodic composite obtained by taking the unit cell boundaries outside = fii 
and almost filling the rest of the unit cell by non-intersecting rescaled and translated 
copies 2 = 2, . . . , n of fi, as illustrated in figure [2l The remainder of the unit cell is 
filled by phase 2 with conductivity o"2- The unit cell structure is periodically repeated to 
fill all space. Let o"c(x) {C for composite) denote this effective conductivity, i.e. in the 
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n 





















Figure 2: A period cell containing rescaled copies of Vt. 

unit cell 

(Tc(x) = o-(x/ai + hi) in fii, i = 1, 2, . . . , iV 

= elsewhere outside Uj fij, (5.1) 

where the scaling constants and translation vectors bj (with ai = 1 and bi = 0) 
are determined by the size and position of each copy fij, so that x/a^ + bj is on the 
boundary of VL if and only if x is on the boundary of fij. Let Pn denote the volume 
fraction in the unit cell occupied by the material with conductivity cr2. Let cr* denote 
the (matrix valued) effective conductivity of this composite, which in general depends 
upon the relative positions of the copies VLi within the unit cell. 
We have the classical variational inequality 

eo • o-*eo < (e^ ■ acGc)^ (5-2) 

which holds for any trial electric field satisfying 

V X e,^ = 0, periodic, (e,^) = eg, (5.3) 

where now the volume averages are over the entire unit cell, rather than just ^2. In 
particular, letting e(x) denote the electric field within VL when the special Dirichlet 
boundary conditions (12. 5p are applied, we may take in the unit cell 

e(^(x) = e(x/aj + b^) in fij, z = 1, 2, . . . , 

= eo elsewhere outside (5.4) 

and periodically extend it. Then we get 
(ec ■ CTc^c) = PnO-2eo ■ eo 
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N 



+ (1 -pn) ^(e(x/ai + hi) ■ a{x/ai + bj)e(x/ai + hi))i 



i=l 



= PnO-2eo ■ eo + (1 -p„)(e(x) ■ o-(x)e(x))n 

= Pn(T2eo ■ eo + (1 - Pn)eo ■ cr^eo, (5.5) 

where denotes an average over the region Qi. Combined with the variational in- 
equahty (15 .2^ this imphes the bound 

eo ■ cr^o < PnO"2eo ■ eo + (1 - p„)eo ■ cr^eo < eo • cr^eo, (5.6) 

where we have used the inequahty cr^ > imphed by (I3.13p . Thus we get 

< < ^""^ (5.7) 

This composite has a volume fraction = (1 — Pn)fi of phase 1. Thus any bound 
"from below" on the effective conductivity cr* , applicable to composites having a volume 
fraction /{ of phase 1, immediately translates into bound "from below" on cr^ . Now 
consider what happens as we increase iV, inserting more and more regions fij, while 
leaving undisturbed the regions f2j already in place, so that p„ — )■ as n — )■ oo. We are 
assured that this is possible. Rescaled copies of any shaped region can be packed to fill 



all space: see, for example. Theorem A.l in Benveniste and Milton (2003). Define 



cr* = lim cr*. (5.^ 



n— >oo 



We are assured this limit exists since if we change the geometry in some small volume 
then the effective conductivity (assuming cxi and a2 are strictly positive and finite) is 



perturbed only by a small amount (Zhikov, Kozlov, and Oleinik 1994). We will call cr 



the effective conductivity tensor of an assemblage of rescaled copies of Vt packed to fill 
all space. Then (15.71) implies 

cr* < cr^, (5.9) 



which is essentially the bound of Huet (1990) applied to this assemblage. Assume the 



bound is continuous with respect to f[ at the point f[ = fi, as expected. Then taking 
the limit n — t- oo the "lower bound" on the effective tensor of composites having volume 
fraction /i must also be a lower bound on cr^. 

In particular, the harmonic mean bound cr* > (cr^^)^^I translates into the elementary 
bound (T^ > {a~^)~^I of Nemat-Nasser and Hori, obtained before. Additionally, in our 
two-phase composite, the effective conductivity cr* satisfies the Lurie-Cherkaev-Murat- 
Tartar bound (Lurie and Cherkaev IWWH \WM IMurat and Tartar 19851 ITartar 19851) 

/iTr[(o-* - ^21)-^] < d/{a, - a2) + /s/as, (5.10) 



[which are a generalization of the bounds of Hashin and Shtrikman (1962 )] where d = 2,3 
is the dimensionality of the composite. Since cr^ > cr* > (J2I it follows that 

(cr* - a2l)-^ > (<T^ - cr2l)-\ (5.11) 

and so (15.101) implies the new bound 

/iTr[(o-^ - a2l)-'] < d/{a, ~ ^2) + /aM. (5.12) 
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By multiplying this inequality by cr| and adding (i/iO"2 to both sides we see that it can 
be rewritten in the equivalent form 

hTr[{a^'l - < d/{a,' - a^') -{d- l)ha^. (5.13) 

As (i^/Tr(A) < Tr(A^^) for any positive definite matrix A we also obtain the weaker 
bound 

lTr[(.^)-i] < a,^ - , (5.14) 

which is a particular case of the universal bounds first derived by Nemat-Nasser and Hori 
(11993[ I1995p . see equation (5.4.9) in their 1995 paper, obtained under the assumption 
that f2 is ellipsoidal or parallelpipedic. (which we see is not needed). 

If one is interested in bounds on the volume fraction /i then f l5.12p implies the upper 
bound 

f < lM + rf/(^l-^2) , . 

l/a2 + Tr[(^^-a2l)-T ^^-^^^ 

To obtain lower bounds on /i, we consider the same periodic composite and apply 
the dual variational inequality 

jo • (O^'jo < (j^ ■ ac%) (5-16) 
valid for any trial current field j satisfying 

Letting j(x) denote the current field within i7 when the special Neumann boundary 
conditions (12.71) are applied, we may take in the unit cell 

j^(x) = j(x/ai + hi) in llj, z = 1,2, . . . ,iV 

= jo elsewhere outside Uj fi,, (5.18) 

and periodically extend it. Substituting this trial field in fl5.16p gives the bound 

«)-^< Wl + (0-', (5.19) 
which in the limit n — )■ oo implies 

cr* > o-^, (5.20) 



which is essentially the bound of Huet (1990) applied to the assemblage of rescaled copies 
of f2 packed to fill all space. 

Thus any bound "from above" on the effective conductivity cr* of composites having 
a volume fraction /( immediately translates into a bound "from above" on (p„(T^^I + 
(cr^)~^)~^. Taking the limit n — )■ cxd and assuming continuity of the bound at /( = /i 
the "upper bound" on the effective tensor of composites having volume fraction /i must 
also be an upper bound on . In particular, the other Murat-Tartar-Lurie-Cherkaev 
bound 

/2Tr[(aiI - < d/{a, - a^) - fi/a^, (5.21) 

implies 

/2Tr[(aiI - o-^)-i] < d/ia, - a2) - h/a,. (5.22) 
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Again using the inequahty (i^/Tr(A) < Tr(A ^) for A > 0, we obtain the weaker bound 

iTr(^^) <a,- — ^'"l (5.23) 

d - 0-2) - /i/(Ti 

which is a particular case of the universal bounds derived by Nemat-Nasser and Hori 
(119931 11995p . see equation (5.3.11) in their 1995 paper, obtained under the assumption 
that Q is ellipsoidal or parallelpipedic (which we see is not needed). 
From f l5.22p we directly obtain the volume fraction bound 

- {Tr[(a,I-.^)-i]}-lM' ^'-'^^ 

giving a lower bound on the volume fraction /i = 1 — /2- 

In the asymptotic limit as the volume fraction goes to zero the volume fraction bounds 
f l5.15p and f l5.24p reduce to those of Capdeboscq and Vogelius fl2003[ I2004p . as shown in 



the two dimensional case by Kang, Kim, and Milton (2011). The paper of Kang, Kim 



and Milton also tests the bounds numerically, and their (two-dimensional) results show 
the bound (I5.15P is typically close to the actual volume fraction for a variety of inclusions 
of phase 1 in a matrix of phase 2. Similarly we can expect that the bound f l5.24p will 
be typically close to the actual volume fraction for a variety of inclusions of phase 2 in a 
matrix of phase 1. 



6 Coupled bounds in two-dimensions 

The tensors cr^ and cr^ obviously depend on ai and cr2, i.e. cr^ = (T^{ai,a2) and 
cr^ = cr^(cri, (72). Let us assume we have measurements of these tensors for an additional 
pair of conductivities {ki,k2), (which could be obtained, say from thermal, magnetic 
permeability, or diffusivity measurements) and let k^ and k^ denote these tensors, 

k^ = (T^(fci, ^2), k^ = (T^ih, k^). (6.1) 

We still let cr^ and cr^ denote the tensors associated with the first pair of conductivities 
(o"i,cr2), with (Ti > (J2- From (15.90 and fl5.20p we have the inequalities 

^2! < cr^ < cr* < cr^ < CTiI, 
A;- < k^ < k* < k^ < A;+I, (6.2) 

where k~ = min{/ci,/c2} and k~^ = max{/ci,/c2} and k* is the effective conductivity the 
composite considered in the previous section when ai and (T2 are replaced by ki and /c2. 
(It can easily be checked that these inequalities still hold if k2 > ki.) 

For two dimensional conductivity from duality (IKeller 1964( Dykhne 1970| ) we know 
the functions cr^ = cr^{ai,a2) and cr^ = cr^(o"i,cr2) satisfy 

(T''{a2,a^) = ^lasRl [o-'^(^i, c^2)]~'R±, 

0-^(^2, ai) = ai(T2Rl[cr^(ai, a2)]~'R±, (6.3) 

where 

Rx=(_°j I) (6.4) 
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is the matrix for a 90° rotation. So if we know these tensors for the conductivity pair 
(fci, /C2); we also know them for the conductivity pair {k2,ki). Hence, by making such an 
interchange if necessary, we may assume without loss of generality that ki > k2, i.e. that 
k~^ = ki and k~ = k2- Finally, by interchanging k with a if necessary, we may assume 
without loss of generality that 

(ri/a2 > k,/k2 > 1. (6.5) 

Optimal bounds on all possible matrix pairs (cr*,k*) for composites having a pre- 
scribed volume fraction /i of phase 1 have been derived by Cherkaev and Gibiansky (1992 ), 
and extended to an arbitrary number of effective conductivity function values by Clark and Milton (1933). 
However it seems difficult to extract bounds on fi from these optimal bounds. Instead 
we consider a polycrystal checkerboard with conductivities 

<t(x) = R^(x)<T*R(x), k(x) = R^(x)k*R(x), with R^(x)R(x) = I, (6.6) 

in which the rotation field R(x) is I in the "white squares" and R^ in the "black squares" . 
By a result of Dykhne (1970) this material has effective conductivities (o"^,I, k^,!) where 

k^, = V det k* . 



V det I 



(6.7) 



Now we replace the "white squares" by the limiting composite considered in the previous 
section (with structure much smaller than the size of the squares) and we replace the 
"black squares" by the limiting composite considered in the previous section, rotated by 
90°. The resulting material is an isotropic composite of phases 1 and 2 and so the pair 
(cr*, k^) satisfies the bounds of Milton (1981b), 

u{k^) < a* < v{k^), (6.8) 

which are attained when the composite is an assemblage of doubly coated disks, where 

2f2a,ia!-al) 



and 



v{k^) 

ai(/c*) 
a2{k^) 



0-2 + 



(/2Cri + /1O-2 + 0-l)(cri + (^2) + {cTi - 02Yai{k^) ' 

Vi^2{al - gj) 

(/2O-1 + /l(T2 + 0-2) ((Ti + 0-2) + (cTi - (J2)2a2(fc*) ' 
(fci + k2)[2f2k,{k, - fc2)/(fci - k,) - if 2k, + hk2 + fci)] 

(fci - k2y 

{k, + k2)[2fMkl - k2)/{k, - k2) - (/2fcl + hk2 + k2)] 

(k, - k2Y 



(6.9) 



(6.10) 



Now for any two symmetric matrices A and B with A > B > we have B ^/^AB > 
I, and so dei{B-'^''^ KQ-^'"^) > 1 implying det(A) > det(B). Thus and ([621) imply 

o"2 < o"Af < 0"* < o"_D < k2 < k^ < k^ < ko < ki, (6-11) 

where we define 

(Tat = a/ det cttv, cxd = a/ det an, k^ = a/ det k^r, kn = a/ det kz). (6.12) 
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The Hashin-Shtrikman bounds ( IHashin and Shtrikman 1962t IHashin 1970p , 



ki 



>k^>k2 + 



(6.13) 



fiki + fik2 + ki /sfci + fik2 + k2 

imply that both ai{k^) and 0:2 (A;*) are non-negative. Hence the denominators in (16. 9p 
are positive and so (16. 8 p imphes 

(cTi - Cr*)[(/2(Ti + /i(T2 + Cri)(0-1 + 0-2) + (CTI - (y2fai{K)] > 2/2(71 ((Ti - ^1), 
(cr* - (T2)[(/2(Ti + /i(T2 + Cr2)(cri + ^2) + (0-1 - a2)^a2(A;*)] > 2/i(T2((Ti - (T2) 

Since aiiko) > ai{k^) and 0:2 (^at) > 0:2 (A;*), we get using (16. lip . 

((Tl - 0-Ar)[(/2Cri + /1O-2 + 0-l)(cri + 0-2) + (o"! - 0-2)^ai (/c^)] > 2/2(71 ((Ji 



(6.14) 



[an - Cr2)[(/20-l + /l(T2 + 0-2) (di + (T2) + ((Tl - (T2)^a2(/i;Ar)] > 2fiCr2{crf 



(6.15) 

As q;i(A;d) and a;2(A;Af) depend hnearly on /i and /2 = 1 — /i, the equations (I6.15P 
readily yield bounds on the volume fraction. Eunjoo Kim has used an integral equation 



solver [as described by Kang, Kim, and Milton (2011)] to compare the bounds (I6.15P 
with the bounds (I5.15P and (I5.24p . Her results are presented in figures [31 HI and 
[51 More numerical results testing the bounds (15.151) and (I5.24p are in the paper by 



Kang, Kim, and Milton (2011). 





Figure 3: The first figure shows the circular body Vt containing an ellipse of phase 1 
surrounded by phase 2. The second figure shows the results for the bounds (I5.15P and 
(I5.24P while the third figure shows the results for the bounds (I6.15p . The bounds are 
for increasing (Ti, with (72 = 1 and (for the third figure) the pairs {ai^ki) are taken as 
(1.1, 1.05), (1.2, 1.1), (1.5, 1.2), (2, 1.5), (3, 2), (5, 3), (10, 5) and (20, 10), with ^2 = ^2 = 
1. Here U{(Ti) and L{(Ti) are the upper and lower bounds on the volume fraction, and 
the true volume fraction is f\ = 0.08. Figure supplied courtesy of Eunjoo Kim. 



7 Coupled bounds in three-dimensions 

We can also derive coupled bounds in three dimensions. Let us assume the phases have 
been labeled so that 

aiki > (72/1:2, i.e. (7i/(72 > k2/ki, (7.1) 
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Figure 4: The same as for figure [3] but with the elhptical inclusion moved closer to the 
boundary of Q. Figure supplied courtesy of Eunjoo Kim. 




Figure 5: The same as for figure [3] but with a non-elliptical inclusion of phase 1 in a 
square region Q. The true volume fraction is /i = 0.0673. Figure supplied courtesy of 
Eunjoo Kim. 
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and by interchanging a with k if necessary let us assume 



ai/(T2 > h/k^. (7.2) 



These two inequahties imply 01/02 > 1 as before. We want to use the inequalities (16. 2 p 
to derive bounds on the volume fraction. As in the two-dimensional case the idea is to 
first construct an isotropic polycrystal, where the polycrystal has the conductivities (16 .6^ 
in which the rotation field R(x) is constant within grains which we take to be spheres. 
These spheres fill all space, and the crystal orientation varies randomly from sphere to 
sphere so that the composite has isotropic conductivities (cr*I, k^T). We use the effective 
medium formula (IStroud 1975} Helsing and Helte 1991 ) which gives 



a* = g{(T^), K = giK), (7.3) 

where for any positive definite symmetric 3x3 matrix A, g = g{A) is taken to be the 
unique positive root of 

+ — — + — ^ — — = 0, (7.4) 



Ai + 2g A2 + 2g A3 + 2g 

in which Ai, A2, and A3 are the eigenvalues of A. This effective medium formula is 
realizable (IMilton 19851 lAvellaneda 1987]) in the sense that it corresponds to a limiting 



composite of spherical grains with hierarchical structure (where any pair of grains of 
comparable size are well separated from each other, relative to their diameter). Note 
that the left hand side of side of (I7.4p increases if any of the eigenvalues Aj increase, 
and decreases if g increases. So g{A) must increase if any or all of the eigenvalues of A 
increase. It follows that g(B) > g{A) if B > A > 0. Hence the inequalities (16. 2p imply 

o"2 < cTAf < 0-* < o-D < 0-1, k^ < k^ < K < kr, < k^ , (7.5) 

where now 

o-jv = 5((crjv), aD=g{(TD), k^ = g{\<iN), fci? = fi'(ko). (7.6) 

We next replace the material in each sphere by the appropriately oriented limiting 
composite considered in the previous section (with structure much smaller than the sphere 
diameter) to obtain a two-phase isotropic composite with {a^,k^) as its conductivities. 
Thus a* must satisfy the upper bound of Bergman (I1976II19781) 

3(72 + (O"! - 0-2)7(fcJ 



where 



and the lower bound 



^.^^^ ^ _ /i/2(fei - ^2) 3A;2 
hki + f2k2-K ki-k2 



y , 3/icr2(ai -a2)(ff2 + 2ai) 

* - ' {f2CTi + h<T2 + 2a2)(cr2 + 2ai) + {a, - a^y^*) ' ^ ' ^ 
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where 



This lower bound was first conjectured by Milton (1981b). A proof was proposed by 



Avellaneda, Cherkaev, Lurie, and Milton (1988) which was corrected by Nesi (1991) and 



Zhikov (1991). 



The lower bound (I7.9P is sharp, being attained for two-phase assemblages of doubly 
coated spheres (IMilton 1981bp . The upper bound fl7.7p is attained at 5 values of '^{k^) 
namely when -f{K) = /a, 8/2/2, 8/2, 8 - 8/i/2, and 8-/1 (IMilton 1981aD . 

The Hashin-Shtrikman bound (IHashin and Shtrikman 1962^ . 

{k. - k2)/{ki - k2) > 3f,k2/{f2ki + hk2 + 2A;2), (7.11) 

implies that I3{k^) is non-negative. Hence the denominator in (17.91) is positive and so the 
inequality implies 

(a*-a2)[(/2fXl + /ia2 + 2a2)((T2 + 2(Ti) + (ori-(T2)'/3(A;,)] > 8/i(T2 (di - (T2 ) (a2 + 2ai) . (7.12) 

The Hashin-Shtrikman bounds can also be rewritten in the form 

f2ki + hk2 + 2k- < //'?^'7^'^; ^ /2A:i + fik^ + 2A;+. (7.18) 

Jiki + f2k2 - k^ 

These inequalities imply 7(fc*) lies between /2 and 8 — /i. Hence the denominator in 
(17. 7p is positive and the inequality can be rewritten as 

(/i^i + 72^2 - ^*)(8a2 + (ai - a2)i{k:)) > f^a^ - ^2)^. (7.14) 

When ki > k2 (ESI) implies Pik^) > I3{k^) and 7(A;d) > l{k^)i and hence 

{ao - o-2)[(/20-i + fia2 + 2(T2)(cr2 + 2ai) + {ai - cr2ff3{kN)] > 8/10-2(0-1 - cr2)(o-2 + 2o-i), 
(/lO-i + /20-2 - o-7v)(8cr2 + (0-1 - o-2)7(/i;d)) > /i/2(o"i - 0-2)^- 

(7.15) 

On the other hand when ki < k2 then (17. 5p implies (5{ko) > /3{k^) and 'jik^) > 7(^*), 
and hence 

(ctd - 0-2)[(/2Cri + /lO-2 + 2o-2)(cr2 + 2cri) + ((Ti - 0-2)^/3(/Cd)] > 8/10-2(0-1 - cr2)(o-2 + 2o-i), 

(/lO-i + /20"2 - criv)(8o-2 + (0-1 - o-2)7(/i;iv)) > /i/2(o-i - 0-2)^ 

(7.16) 

Since I3{kj^) and /3{k£)) depend linearly on the volume fractions /i and /2 = 1 — /i, the 
first inequalities in (I7.15P and (I7.16P also depend linearly on the volume fraction and 
easily yield bounds on the volume fraction. On the other hand, finding bounds on the 
volume fraction from the second inequalities in (I7.15P and (17.160 , involves solving a cubic 
equation in /i. So instead of analytically computing the roots of this cubic it is probably 
better to numerically search for the range of values of /i where the second inequalities 
in (I7.15P and (17. 16p are satisfied. 



18 



8 Bounds for elasticity 

Let us consider solutions to the linear elasticity equations 

t(x) = C(x)e(x), V-e = 0, e = (Vu + (Vu)'^)/2, (8.1) 

within fl, where u(x), e(x) and t(x), are the displacement field, strain field, and stress 
field, and C(x) is the fourth order elasticity tensor field 

C(x) = x(x)Cl + (l-x(x))C^ (8.2) 

in which and are the elasticity tensors of the phases, assumed to be isotropic with 
elements, 

^ijke = f^hi^ikSje + 6ii6jk) + {kh - 2^h/d)5ij5ki-, h = 1,2, (8.3) 

in which ci = 2 or 3 is the dimensionality, and /ii,/i2 and ki, K2 are the shear and bulk 
moduli of the two phases. From boundary information on the displacement uo(x) = u(x) 
and traction f(x) = t(x) ■ n we can immediately determine, using integration by parts, 
volume averages such as 

l"l JdQ 

(6) = ^ / (nu^ + un^)/2, 

l"l JdQ 

(r) = ^ / xf^, (8.4) 



in which " : " denotes a contraction of two indices. 

There are two natural sets of boundary conditions. For any symmetric matrix eo we 
could prescribe the special Dirichlet boundary conditions 

u(x) = eox, for x G dil, (8.5) 

and measure tq = (r). Here, according to (18. 4p . eo equals (e). Since tq is linearly related 
to eo we can write 

To = C^eo, (8.6) 

which defines the elasticity tensor cr^ {D for Dirichlet). Alternatively for any symmetric 
matrix To we could prescribe the special Neumann boundary conditions 

r(x) ■ n = To • n, for x G dVl, (8.7) 

and measure eo = (e). Here, according to (18. 4p . tq equals (r). Since eo is linearly related 
to To we can write 

eo = (C^)-Vo, (8.8) 

which defines the elasticity tensor cr^ (D for Dirichlet). It is easy to check that and 
satisfy all the usual symmetries of elasticity tensors. 
Directly analogous to (I3.13P we have the bounds 

(C-^)-^ < (C), {C-y^ < C"" < (C) (8.9) 
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of (INemat-Nasser and Hori 1993p . and directly analogous to fl3.14p for any boundary 
condition (not just the special boundary conditions flS.Sp and (18. 7p ) we have the bounds 

(e ■ r) > eo ■ a^eo, (e ■ r) > To(o-^)-Vo, (8.10) 



where eo = (e) and To = (r), due to Willis in 1989 private communication to Nemat- 



Nasser and Hori and presented by Nemat-Nasser and Hori (1993). 
Also directly analogous to (15. 9 p and (I5.20p we have the bounds 

<e* < c^, 



(8.11) 



where C* is the effective elasticity tensor of any assemblage of rescaled copies of fl packed 



to fill all space. These are essentially the bounds of Huet (1990) applied to this as- 
semblage. Thus "lower bounds" on C* directly give "lower bounds" on and "upper 
bounds" on C* directly give "upper bounds" on C^. In particular, in two dimensions lower 



and upper bounds on eo-C*eo have been obtained by Gibiansky and Cherkaev (1984) (for 



the equivalent plate equation) and also by Allaire and Kohn (1993a). Assuming that the 



phases have been labeled so that /ii > /i2 (/^i — ^2 could be either positive or negative) 
and letting ei and €2 denote the two eigenvalues cq, the bounds imply 



eo ■ C^eo > (ei + e2)7(/i/'^i + 12/^^2) + (ei 

if |ki - /t2|(/l/i2 + /2/il)|ei + 62! < 



e2)V(/l//il + /2//i2), 

/i2|(/ifi:2 + /2ft:i)|ei - 62! 



eo ■ C^eo > (ei + e2)l/i/ti + /2/t2) + (ei - £2)1/1/^1 + /2/U2) 
_ J. J. [1^1 - /t2||ei + £2! + 1/^1 - /^2||ei - £2!]^ 

if (/i2 + /1K2 + /2/ti)|ei - £2! > /2|/ti - fi;2||ei + €21 

and |ki - K2|(/i/U2 + /2/ii)|ei + £2! > - fi2\{fii^2 + f2Hi)\€i - €2! 



and 



eo ■ C^eo > I22{ei - £2)' + 



2 K1K2 + At2(/l/tl + /2/«2) 



ei + 62 



1^2 + fll^2 + /2/«l 

if (/i2 + /i/t2 + /2/ti)|ei - £2! < f2Wi - /«2||ei + £2!; 

eo • C^eo < (ei + e2)'(/i/€i + f2^^2) + (d - e2)'(/i/^i + /2/i2) 
_ [Wi - ^2||ei + ^2! - l/ii - /i2||ei - 62!]^ 

/l(/i2 + ^2) + /2(/il + fi^l) 

if - K2||ei + £2! < (/Ui + /1K2 + /2'ti)|ei - €2! 

and - /i2||ei - £2! < (/t+ + /1/U2 + /2/ii)|ei + 62!; 



(8.12) 



eo-C eo < Aii(ei - 62) H — . — (ci + £2) , 

yUi + /i/^2 + /2«:i 

if fi\f^i - «;2||ei + £2! > (/ii + /i«2 + /2«;i)|ei - £2!; 



r^N ^ / I n2 I /^l/^2 + '«+(/lA*l + 72/^2) / n2 

eo ■ <^ eo S + '^2 H —7 —7 l^i — €2 , 

/«+ + /i/i2 + /2yUi 

if f+\ni - /i2||ei - £2! > (/t+ + /1/U2 + /2/ii)|ei + £2!, 



.13) 
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where k+ is the maximum of ki and K2 and /+ is the volume fraction of the material 
corresponding to k,^. 

The corresponding three-dimensional bounds follow directly from ( IS.lip and the 
bounds of Allaire and Kohn (1993b), but are not so explicit. Assuming that the Lame 
moduli 

Ai = Ki — 2yUi/3 and X2 = ^2 — 2/i2/3 (8-14) 

of both phases are positive, and that the bulk and shear moduli of the two phases are 
well-ordered with 

Ki > K2 > and /ii > yU2 > 0, (8.15) 

these bounds are 

eo : C^eo > eo : C2eo + /i max[2eo : 77 - 77 : (Ci - C2)"^77 - /25'(^)], 

60 : C^eo > eo : CiEq + /2 min[2eo : f] + f] : (Ci - C2)-''n - fiHr))], (8.16) 

where (7(77) and h{r}) are function of the eigenvalues rii,ri2, and rj^ of the symmetric matrix 
T]. Assuming that these are labeled with 



r]i<r]2< r]3, (8.17) 



we have 



and 



4/i2 4(A2 + /i2) 2(A2 + /i2) 

/X .f ^ A2 + 2/i2 . , . 

A2 + 2;U2 2(A2 + yU2) 

Kv) = inm{ril rjl rjl}. (8.19) 

The bounds f l8.12p . f l8.13p and fl8.16p can be used in an inverse way to bound the 
volume fraction fi = 1 — /2, for a single experiment when for special Dirichlet conditions 
eo is prescribed and To (= C^^o) is measured, or when for special Neumann conditions 
To is prescribed and eo (= C^eo) is measured. Allaire and Kohn (1993b) also derive 
bounds on the complementary energy and these imply 

To : (C^)-Vo > To : CrVo + /2max[2To : C'C : (C,' -C^Y^C- fiC : e^C + fMC.C)], 

(8.20) 

and 

To : (O^Vo < To : C2-Vo + /imm[2To : C + C : {C2' -C^T'C- f2C ■ C2C + f29{C2C)]- 

(8.21) 

The bound in f l8.20p is particularly useful when To = —pi, corresponding to immersing 
the body f2 in a fluid with pressure p. Then from measurements of the resulting volume 
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change of the body one can determine tq : (C^)^^ro = — pTreo. Let us assume Ai > 
and set <^ = al + A, with A being a trace free matrix with eigenvalues ai, 02 and 03. 
Then we have t] = CiC = 2/ii(A;I + A) where k = a[l + 3Ai/(2/ii)]. Substitution gives 

[C:CiC-/^(CiC)]-«'[I:CiI-/i(CiI)] 



2/ii 



222 
+ + CI3 



2^1 



Ai + 2/ii 

> 2/xi[a^ + al + al — min{2aiA; + a^, 2a2k + a\, 2a3k + Og}] 



min{(A; + aif - k^, {k + 02)^ - A;^ {k + a^f - k^} 

(8.22 



which is surely positive since min{2aifc + ai, 2a2k + a^, 2a^k + a\} < a| where j is such 
that kttj is non-positive. (Note that ai, 02 and ag cannot all have the same sign since they 
sum to zero). Consequently when tq = —pi the maximum over C, in fl8.20p is achieved 
when A = and taking the maximum over a gives 



— p Tr eo > 



Ki 



+ 



/2 



+ 



Ki— 3/ti+4/ii 



(8.23) 



or equivalently 



- pI (Tr eo) < k^^^ = 1^1-^ 



/2 



/ {1^1 - H2) - fi/ (f^i + 4/ii/3) ' 



(8.24) 



where k]jsh upper bulk modulus bound of Hashin and Shtrikman (1963) and 



Hill (1963). The inequality f l8.23p can be rewritten as 



/2 



< 



K.1K2 



-{Tt eq/p) - 1/ki ~ Ki - K2 ' 3Ki + 4/ii 



+ 



(8.25) 



where we have used the fact that — (Treo/p) — l/^i is positive (since Tq : (C ) Tq > 
Tq : (Ci)~^To by (18. 9p ). This then yields the volume fraction bound 



/2< 



+ 



3ki+4/^i 



1 



-(Treo/p)-l/Ki ' 3Ki+4Ati 



+ 



(8.26) 



which we expect to be closest to the actual volume fraction when phase 2 (the softer 
phase) is the inclusion phase. Thus the bound may be particularly effective for estimating 
the volume of cavities in a body. Note that if some granules of phase 1 lie within these 
cavities, then such granules will not contribute to this volume fraction estimate, but will 
contribute to the overall weight. If the weight of the body has been measured (and the 
density of phase 1 is known) this provides a way of estimating the volume of granules of 
phase 1 which lie within the cavities. 

When multiple experiments have been done, and the full tensor or has been 



determined, then the "trace bounds" of Zhikov flT988l [T99T]) and [Milton and Kohn (1988[ ) 
can be used. (These generalize the well known Hashin-Shtrikman (I1963p bounds to 
anisotropic elastic composites.) Define the two traces 



Tr,, A = Aiijj/d, TtsA = A 



i^ajj/ d), 



(8.27) 
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for any fourth order tensor A. with elements Aijke in spatial dimension d. Then, assuming 
the moduli of the two-phases are well ordered satisfying f lS.lSp . their lower and upper 
"bulk modulus type bounds" imply, through (IS.lip . the universal bounds 



/iTr,[(C 



D 



< 



1 



/2 



f2TTh[{ei-c 



N\-l-\ 



< 



d(Ki — K2) dK2 + 2[d — l)yU2 
1 /l 



f8.28) 



d{Ki — K2) dKi + 2{d — l)/ii ' 
while their lower and upper "shear modulus type bounds," imply the universal bounds 

{d-l){d + 2) d{d-l){K2 + 2fi2)f2 



hTT,[{C^-C2)-'] < 



/2Tr,[(Ci-C 



< 



4(/ii - /i2) 2n2{dK2 + 2{d - 1)112) 
(d-l)(rf + 2) + 2/ii)/i 



4(/ii-/i2) 2/ii(rfKi + 2(rf- 



^.29) 



Since these inequalities depend hnearly on /i = 1 — /2 they can easily be inverted to 
obtain bounds on fi given or C^. 

As noted by Milton and Kohn (1988 ) the lower and upper "bulk modulus type bounds" 



are tighter than those obtained by Kantor and Bergman (1984) and Francfort and Murat (1986), 
which imply 



Tih < dKi 



f2 



h 



l/Tr,[(C^)-i] > dK2 + 



fl 



1 



+ 



/2 



^.30) 



K2) dK2+2{d—l)iM2 



For bodies fl of ellipsoidal or parallelopipedic shape the universal bounds fl8.30p were 
obtained by Nemat-Nasser and Hori ( I1993[ [T995|) : see the equations (4.3.9) and (4.4.8), 
with / = 1, in their 1995 paper. Their other bounds, with 1 = 2, which incorporate the 
"shear responses" of the tensors and are improved upon by the bounds ( 18.29^ as 
can be seen using the inequality 



TtsA-' > (rf-l)2(d + 2)V(4Tr,A), 
which holds for any positive definite fourth-order tensor A. 



(8.31) 
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